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Abstract 

We present some simple methods to find gluon distribution from 
analysis of deuteron F% structure function data at moderately low- 
x. Here we use the leading order(LO) Altarelli -Parisi(AP) evolution 
equation and New Muon Collaboration (NMC) deuteron F2 structure 
function data to extract gluon distribution. We also compare our re- 
sults with those of other authors. 
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1. Introduction: 

The measurements of the proton and the deuteron F 2 (x, Q 2 ) structure 
functions by Deep Inelastic Scattering (DIS) processes in the low-x region 
where x is the Bjorken variable have opened a new era in parton density 
measurements [1]. It is important for understanding the inner structure of 
hadrons as well as examining QCD, the underlying dynamics of strong inter- 
action among the partons inside hadrons. It is also important to know gluon 
distribution G(x, Q 2 ) inside hadron at low-x because gluons are expected to 
be dominant in this region. On the otherhand, measurement of gluon distri- 
bution directly from experiments is a difficult task. It is, therefore, important 
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to measure gluon distribution G(x, Q 2 ) indirectly from proton or deuteron 
structure functions F 2 (x, Q 2 ). Here the representation for the gluon distribu- 
tion G(x) = xg(x) is used where g(x) is the gluon density. A few number of 
papers have already been published [2-9] in this connection. Here we present 
two alternative methods to relate Gluon distribution G(x, Q 2 ) with deuteron 
F 2 (x, Q 2 ) structure function and their differential coefficients with respect to 
InQ 2 and x, i.e. dF 2 (x, Q 2 )/dlnQ 2 and dF 2 (x, Q 2 )/dx for fixed values of Q 2 . 
We report for the first time some methods to extract gluon distributions from 
deuteron F 2 structure function data. Our methods are simpler with less ap- 
proximation and more transperent. Of course, there exist some established 
methods [10] for extracting gluon distributions from data based on global 
fits. In these methods, momentum distribution and other constraints [11] 
are used to get gluon distributions. But our methods are based on the direct 
solutions of QCD evolution equations which may be some good alternatives. 

Section 1 in our paper is the introduction. Section 2 deals with the theory 
for extracting gluon distribution from deuteron F 2 structure function data. 
Section 3 is the result and discussion and section 4 is the summary and con- 
clusion. 

2. Theory: 

In the LO analysis deuteron structure function is directly related to singlet 
structure function [12]. On the otherhand, the differential coefficient of sin- 
glet structure function F 2 s (x ) Q 2 ) with respect to InQ 2 , i.e. dF^x, Q 2 )/dlnQ 2 
has a relation with singlet structure function itself as well as gluon distribu- 
tion function [12] from AP evolution equation [13-16]. Thus it is possible 
to calculate gluon distribution from singlet structure function or ultimately 
from deuteron structure function also. 

The LO AP evolution equation for singlet structure function [12] is given 



by 
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where, t = ln(Q 2 /A 2 ) and A f = 4/(33 - 27V/), N f being the number of 
flavours and A is the QCD cut off parameter. Now, 



oo 

X 



.E^s+sj:** (2) 

1 Z fc=0 fc=l 

We have, 1 — :e<2:<0=^|,2:|<1 which implies that the expansion (2) is 
convergent. Now by the Taylor expansion [17] we get, 

n (^ z ,^ F - ix , t)+x p a Ji£A (3) 

andC?(^,i)^M) + zf>^|^ (4) 

neglecting the higher order terms. 

But as a matter of fact, we can not neglect the higher order terms as these 
terms are not small in Regge-like [7,18] or in Double-logarithmical [7,19] be- 
haviours for singlet structure function or gluon distribution function. On 
the otherhand, it has been shown that this Taylor expansion method is suc- 
cessfully applied in calculating (Revolution [20-22] or re-evolution [23] of 
structure function with excellent phenomenological success. Some authors 
[3-5] again applied this method to extract gluon distribution from proton 
structure function. It was suggested that [23] one possible reason for success 
of this method may be due to simplification of QCD processes at low-x for 
momentum constraints. 

Putting equations (3) and (4) in equation (1) and performing z- integrations 
we get 

2 A. - -/-[A s (x)Fi(x, t) + B s (x)G(x, t) + C s (x) 2 v ' 
at t ox 

+flsW ^) ] = , (5) 

where, A s (x) = 3 + 4Zn(l - x) + 2{(1 - x)(-2 + (1 - x)/2)}, 
B s (x) = (3/2)iV / {(l -x)(x+ (2/3)(l - x) 2 )}, 
C s (x) = 2x{ln(l/x) + (1 - x)(l - (1 - x)/2)} 
and D s {x) = (3/2)N f {ln(l/x) - (1 - x)(l + (2/3)(l - x) 2 )}. 

Now, we can apply two methods to extract gluon distributions: 
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First method: 

At very low-x limit, x — > 0, the functions A s (x), C s (x), and D s (x) become 
vanished and B s (x) = Nf. Equation (5) then becomes simplified and we get 

^G(xt) = -L*^ (6 ) 
AfN/ dt ■ [) 

Equation (6) is a very simple relation between gluon distribution function 
with the differential coefficient of singlet structure function with respect to 
t. 

Second method: 

Recasting equation (5) we get 

C( t\ i D s( x )dG(x,t) _ 1 dFj(x,t) 
{X,) + B s (x) dx A f B s {x) dt 

A s (x) C s (x)dFj(x,t) 



B s (x) 2X ' ' B s (x) dx 

Now D s {x)/ B s {x) is very small at low-x, lim x ^ D s (x) / B s (x) = 0. So, apply- 
ing the Taylor expansion series we can write 

Thus equation (7) gives 

G(x',t) = Kl (x)X.?^^ + K 2 (x)?^^ + K 3 (x)^(x,t), (8) 

, D.Jx) T ^ . . 1 

where, x = x + , Ki(x) = 



B s {xY iV ' A f B s 



B s (x) ° v ' B s (xY 

Equation (8) is also a simple relation between gluon distribution function 
with the differential coefficients of singlet structure function with respect to 
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t and x, and with singlet structure function itself. If we try to combine the 
last two terms of equation (8) let us take common K 3 (x) from both the terms 
which reduce to 



K 3 (x) 



F°(x,t) + 



K 2 (x) dFj(x,t) 



KM 



dx 



But K 2 (x)/K s (x) 1S n °t small at low-x and therefore these two terms can 
not be combined to one as in the case of gluon by applying Taylor expansion 
series. 

The relation between deuteron and singlet structure function at LO [12] 

is 



Then we get 



F*(x,t) = ^F 2 s (x,t) F°(x,t) = ^F 2 d (x,t). 



dF£(x,t) _ 9 dF$(x,t) 
~ 5' 



and 



dt 5 dt 
dF$(x,t) _ 9 dF$(x,t) 
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dx 



dx 



(9) 
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Putting equations (9), (10) and (11) in equations (6) and (8), we get repec- 
tively 



G{x,t) 



9t dF$(x,t) 



5A f N 



/ iV / 



dt 



G(x',t) = 



, « dF$(x,t) .dF$(x,t) v„ rf/ , 

KM-t- I + KM Q + K 3 (x)F^(x, t) 



dt 



dx 



and (12) 
(13) 



which are our main results. From these equations it is seen that if we have 
deuteron structure function and their differential coefficients with respect to 
t and x at any x for a fixed value of t = t , we can calculate gluon distribution 
function at x (first method) from equation (12) or at x' = x + D s (x)/B s (x) 
(second method) from equation (13) as a LO analysis. 

For analysis of our results, we use NMC 15-parameter function [24,25] 
which parametrized their data for proton and deuteron structure functions 
for Q 2 values from 0.5 GeV 2 to 75 GeV 2 and low-x values from 0.002 to 0.6. 
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This parametrization can also well describe the SLAC [26] and BCDMS [27] 
data, and Fermilab [28] low-x data. The function used to describe proton as 
well as deuteron data is given by 



F 2 {x,Q 2 ) = A{x) 



ln(Q 2 /A 
[HQl/A*)_ 



1 + 



C(x) 

Q 2 



(14) 



Here Q 2 = 20 GeV 2 and A = 0.250 GeV, 



A(x) 



x ai (l 



{a 3 + a 4 (l -x) + a 5 (l - x) 2 + a 6 (l - x) 3 + a 7 (l - x) 4 }, 



B(x) — bi + b 2 x + b 3 /(x + 6 4 ) and 
C(x) = c\x + C2X 2 + C3X 3 + C4X 4 

where a±, a 2 , 0,3, a 4 , a$, a^, 0,7, b±, b 2 , 63, 64, c±, c 2 , C3 and C4 are the 
15-parameters used to fit the data. Actually two different sets of these pa- 
rameters are used to describe proton and deuteron structure functions in the 
same equation, equation (14). Thus for the respective sets of parameters 
equation (14) gives the deuteron structure function as 



F*(x,t)=A(x) 



- 1 - 


B(x) 




To. 




1 + £.C(x) 



(15) 



where t = ln(Q 2 /A 2 ) and t = ln(Ql/A 2 ). Differentiating F$(x,t) with 
respect to t and x we get respectively 



dFf(x,t) (B{x) 



dt 



dF 2 d (x,t) 
dx 



\ 1 
dA(x) 
dx 



l)F*{x,t)+A{x) 



,t\B(x) 
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+ A(x) 
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and 

dB{x) 
dx 
e~ l dC(x 



(16) 



dx 



(17) 



where, 



dA(x) 
dx 



x 



a 2 



A(x) -x ai (l-x) a2 .{a 4 + 2a 5 (l-x) 



1-x, 

+3a 6 (l - x) 2 + 407(1 - x) 3 }, 
dB(x) 



dx 



b 2 - 



[x + b 4 y 



and 
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- = ci + 2c 2 x + 3c 3 x 2 + 4c 4 x 3 . 
ox 

Now putting equations (15), (16) and (17) in equatins (12) and (13) we 
can easily calculate gluon disatributions at x (first method) or x' = x + 
D s (x)/B s (x) (second method) respectively. 



3. Results and Discussion: 

The NMC 15-parameter function [24,25] we use, parametrized the NMC 
for Q 2 values from 0.5 GeV 2 to 75 GeV 2 and low-x values from 0.002 to 0.6 
which also well describe the SLAC [26] and BCDMS [27] data and Fermilab 
[28] low-x data. As the data range of x we use is moderatelt low, we will 
restrict our analysis for Q 2 values from 10 GeV 2 to 60 GeV 2 and low-x values 
from .1 to 0.001. We can not extend our analysis to HERA low-x region [1] 
due to lack of deuteron F 2 structure function data in that region. 

In Fig. 1(a) and Fig. 1(b) gluon distributions obtained by our first method 
(equation (12)) from NMC deuteron parametrization from the 15-parameter 
function are presented at Q 2 = 10 GeV 2 and 60 GeV 2 respectively. The mid- 
dle lines are the results without considering the error. The upper and the 
lower lines are the results with parameter values by adding and subtracting 
the statistical and the systematic errors with the middle values respectively. 
It has been seen that the middle lines almost coincide with the upper ones. 
We calculate gluon distributions for x-values from 10 _1 to 10~ 3 for both 
Q 2 = 10 GeV 2 and Q 2 = 60 GeV 2 . In both the cases G(x, Q 2 ) values in- 
creases when x decreases as expected, but G(x, Q 2 ) is higher in Q 2 = 60 GeV 2 
than in Q 2 = 10 GeV 2 for same x especially in lower-x side. Moreover, rate 
of increment of G(x, Q 2 ) is very high from x = 10" 1 to 10~ 2 . But the rate 
decreases to some extent to lower-x region. 

Exactly in the similar way, in Fig. 2 (a) and Fig. 2(b) gluon distributions 
obtained by our second method (equation (13)) from NMC deuteron parame- 
trization from the 15-parameter function are presented at Q 2 = 10 GeV 2 
and 60 GeV 2 respectively. All discussion are exactly same as for Fig. 1(a) 
and Fig. 1(b). But overall values of G(x, Q 2 ) are higher in second method 
than in first one for any value of x. For example, G(x, Q 2 ) medium values 
are almost 20% and 25% higher in second method than in first method for 
Q 2 = 10 GeV 2 and Q 2 = 60 GeV 2 respectively at x — 10~ 3 . This is because 
in our first method we apply very low-x approximation and neglected A s (x), 
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C s (x) and D s (x) in equation (5) as they are vanishingly small at very low-x 
to obtain equation (6) and then equation (12). On the otherhand, in our 
second method we do not apply such approximation and automatically the 
contributions from these functions have been included in equation (13). 

In Fig. 3, comparison of gluon distributions obtained by Sarma and Medhi 
method (SM), Bora and Choudhury method (BC), Prytz method, our first 
method (SA 1st) and our second method (SA 2nd) is presented for middle val- 
ues only for Q 2 = 60 GeV 2 . Values are higher for the results of other authors 
with proton structure function data than of ours with deuteron structure 
functions data. This is actually due to the fact that the scaling violations of 
deuteron structure functions F$(x, Q 2 ) with respect to lnQ 2 (= t) are them- 
selves considerably less than those of HERA proton data due to HI [29,31] 
and ZEUS [30,32] collaborations and these scaling violations are directly 
proportional to gluon distributions in the formulas used by BC and Prytz to 
calculate gluon distributions. These HERA proton data covers x values upto 
at least ~ 10~ 4 in comparison with those of NMC data which covers upto 
~ 10~ 3 only. 

Gluon distributions increase as x decreases due to all the authors as ex- 
pected from QCD analysis. Moreover, gluon distribution by our first method 
is lowest and by Sarma and Medhi method is the highest for a particular 
low-x. 

4. Summary and Conclusion: 

In this article we present for the first time a method to extract gluon dis- 
tribution from the measurement of moderately low-x deuteron structure func- 
tions and their differential coefficients with respect to t = InQ 2 and x. Here 
we use LO AP evolution equation to relate gluon distribution function with 
moderately low-x structure functions or their differential coefficients. In our 
analysis we use only NMC deuteron data parametrization by a 15-parameter 
function. We find that gluon distribution from deuteron also increases when 
x decreases as in the case of proton as usual. We can not compare our re- 
sult of NMC data with other because low-x deuteron data is not sufficiently 
available. Moreover, no other work to calculate gluon distribution function 
from deuteron data has been so far reported. But we compare our result with 
gluon distributions due to other authors Sarma and Medhi, Bora and Choud- 
hury and Prytz calculated from low-x proton data. We see that our result is 
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to some extent less as differential coefficient of deuteron structure function 
with respectr to t = InQ 2 is much less than that of proton structure function. 

In our method the first order approximation in Taylor expansion in — 
z),t) and G(x/(1 — z),t) is used, i.e. only terms having first order differen- 
tial coefficients dF^x, t)/dx and dG(x, t)/dx are used. Scope is still there to 
include higher order terms of the Taylor expansion series. Of course, our pre- 
liminary work including second order differential coefficients (9 2 F 2 s (x, t)/dx 2 
and d 2 G(x, t)/dx 2 reveals that they could not contribute significantly. More- 
over, this is only a LO analysis. We can expect better result if we include 
next-to-leading order (NLO) and subsequent terms in perturbative QCD. 
Work is going on in this regards. 
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Fig. 1(a) and Fig.l(b):Gluon distributions obtained by our first method 
(equation (12)) from NMC deuteron parametrization from the 15-parameter 
function are presented at Q 2 = 10 GeV 2 and 60 GeV 2 respectively. The mid- 
dle lines are the results without considering the error. The upper and the 
lower lines are the results with parameter values by adding and subtracting 
the statistical and the systematic errors with the middle values respectively. 
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Fig. 2 (a) and Fig. 2(b): Same result as in Fig. 1(a) and Fig. 1(b) by our sec- 
ond method (equation (13)). 
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Fig.3:Comparison of gluon distributions obtained by Sarma and Medhi method 
(SM), Bora and Choudhury method (BC), Prytz method, our first method 
(SA 1st) and our second method (SA 2nd) for middle values only without 
considering any error. 
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